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Abstract 

The simultaneous presence of seemingly incompatible properties of solids often provides a unique 
opportunity to address questions of fundamental and practical importance. The coexistence of 
ferroelectric and topological orders is one such example. Ferroelectrics, which have a spontaneous 
macroscopic polarization switchable by an applied electric field, usually are semiconductors with 
a well-developed wide band gap mm with a few exceptions [IHS]. On the other hand, time- 
reversal symmetric Z 2 topological insulators (TI) [7], characterized by robust metallic surface states 
protected by the topology of the bulk, usually are narrow-gap semiconductors (< 0.7 eV) [HI [9] 
which allow band inversion induced by the spin-orbit interaction. To date, a ferroelectric topological 
insulator (FETI) has remained elusive, owing to the seemingly contradictory characters of the 
ferroelectric and topological orders. Here, we report that the FETI can be realized in halide 
perovskite CsPbIs under strain. Our first-principles study reveals that a non-centrosymmetric 
ferroelectric structure of CsPbIs is energetically favored under a wide range of pressures, while 
maintaining its topological order. The proposed FETI is characterized by switchable polar surfaces 
with spin-momentum locked Dirac cones, which allows for electric-field control of topological surface 
states (TSSs) and the surface spin current. Our demonstration of a FETI in a feasible material 
opens doors for future studies combining ferroelectric and topological orders, and offers a new 
paradigm for diverse applications in electronics, spintronics, and quantum information. 
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A TI is a topological state of matter that is characterized by metallic boundary modes 
(edge or surface) arising from the non-trivial topology of an insulating bulk electronic struc¬ 
ture [7] . The metallic nature of the boundary states is protected by time-reversal symmetry, 
providing transport channels that are robust against ordinary non-magnetic impurities. Fun¬ 
damental scientihc interest and potential applications of TIs in electronics, spintronics, and 
quantum information have driven tremendous efforts to hnd robust material realizations of 
TIs. To date, more than 50 candidate materials ng have been suggested theoretically in 
diverse material classes with some confirmed experimentally, including chalcogenides ffu- 
[TH], oxides 0 HBHIH], ternary Heusler compounds [E], halide perovskites [201 El], and a 
non-centrosymmetric Rashba semiconductor [22] . A characteristic feature common to these 
materials is a narrow band gap originating from the band inversion mechanism. Without 
the spin-orbit interaction, these materials would be normal insulators (NIs); the spin-orbit 
interaction interchanges the characters of the valence and conduction bands, resulting in a 
non-trivial topology. For this to occur, the band gap must be smaller than the strength of 
the spin-orbit interaction, typically less than 1 eV, in order to host TIs. On the other hand, 
ferroelectricity is usually found in transition metal oxides with wide band gaps, due to the 
large difference in electronegativity between oxygen and transition metals. Most technolog¬ 
ically important ferroelectric materials (e.g., PbTiOs and BaTiOa) have band gaps larger 
than 3 eV [2]. The lowest band gap for a ferroelectric oxide is 1.39 eV, obtained recently in 
the single-phase solid oxide solution [KNb03]i_a;[BaNii/2Nbi/203_5]a;(KBNN0) [ 1 ]. 

While the observed narrow band gap in various TIs seems to exclude the co-existence of 
ferroelectricity, there is no fundamental physics preventing the realization of both ferroelec¬ 
tric and topological order in a single-phase material. A ferroelectric topological insulator 
that combines switchable polarization and robust TSSs is likely to support novel device 
functionalities. Here, we demonstrate through first-principles calculations that the FETI 
can be realized in cubic halide perovskite CsPbIs under strain in the absence of Pble octa¬ 
hedron rotation. The consequences of the coexistence of ferroelectricity and the non-trivial 
topology are explored. We find that the TSSs in FETI are tunable by the electrostatic 
potential associated with the polarity of the surfaces, which allows for simplified nanoscale 
engineering of topological p-n junctions [2S] when combined with domain engineering. [21] 

Using density functional theory (DET) calculations (see Methods), we find that the cubic 
halide perovskite CsPbIs experiences a paraelectric-ferroelectric structural phase transition 
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under strain. In the absence of Pble octahedron rotation, the cubic CsPbIs with ferroelectric 
distortions (Pb atom displacing away from the center of Pble octahedron) has lower forma¬ 
tion energy (~ 10 meV per unit cell) than the paraelectric CsPbIs for a wide range of lattice 
constants (a < 6.05 A) at 0 K. Free energy calculations including zero-point energy and 
the entropy contribution from phonons at room temperature also indicate that the system 
should undergo a ferroelectric distortion when the lattice constant o is compressed smaller 
than ~ 5.90 A. 

Pressure also induces a topological phase transition from a normal insulator to a topologi¬ 
cal insulator in CsPbIs. Yang et ah and Jin et al. predicted and have shown that paraelectric 
CsPbIs is a topological insulator under hydrostatic pressure [201 El]. We find that the elec¬ 
tronic structure of ferroelectric CsPbIs at a < 6.04 A is adiabatically connected to the 
topological phase of the paraelectric CsPbIs. More specifically, as shown in Fig. Ic, the 
band gap remains open during a structural deformation from the paraelectric topological 
phase to the ferroelectric phase at a — 5.79 A. The bulk band structures of ferroelectric 
CsPbIs, calculated with varying a, also show evidence of the topological phase transition 
via gap closure. The bulk gap closes at a = 6.04 A(Fig. lb) , exhibiting a nodal point 
along the R-X line, and reopens at a < 6.04 A. Interestingly, a nodal point occurs not 
only along the high-symmetric direction, but also along an arbitrary direction away from 
R on the kz = tt/o plane, forming a line node encircling R. As previously discussed [25] . 
it is a characteristic feature of topological phase transition between NI and TI in three- 
dimensional non-centrosymmetric systems that the gapless states, in general, appear off 
the time-reversal invariant point (R) in one-dimensional parameter subspace of the (k, m) 
space, where m is an external parameter tuning a topological phase transition. Whereas 
the specific geometry of nodal points can depend on the details of crystalline symmetries, 
their formation again supports a ferroelectric topological phase at a j 6.04 A. We directly 
calculate the Z 2 topological invariants m, employing the method proposed in Ref. EHl The 
indices are (/ro!/^i/^ 2 /^ 3 ) = (RIH), further confirming that ferroelectric CsPbR is a strong 
topological insulator under high pressure. As band gaps are systematically underestimated 
within DFT, electronic correlations introduced by higher level theories such as the GW ap¬ 
proximation are expected to increase the band gap and hence lower the critical value of a at 
which CsPbIs becomes a TI[2l|. Since strain favors the ferroelectric phase (in the absence of 
Pble cage rotation), CsPbR is likely to be a FETI even in the presence of electron-electron 
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correlations. 

To illustrate the TSSs of a FETI, we have calculated the surface band structures using 
a slab model (see Methods). Ferroelectric materials possess three different types of sur¬ 
faces depending on the orientation of the polarization: non-polar positive-polar and 
negative-polar P~ surfaces. We first investigate non-polar surface states by exposing the 
(001) surfaces with in-plane [100] polarization (Fig. 2a). As shown in Fig. 2b, the surface 
states exhibit a Dirac-cone-like energy dispersion at the M point (0.5,0.5) of the surface 
Brillouin zone, confirming the nature of TSSs. 

The surface states of the polar surfaces are studied using a (100) slab with [lOOj-oriented 
polarization (Fig. 3a). Figure 3b shows the surface band structure, with states from 
(red) and P~ (blue) surfaces resolved separately for bands close to the Fermi energy (see 
Methods). The out-of-plane polarization breaks the degeneracy of the surface states, shifting 
upward in energy the Dirac cone (p-type) associated with the P~ surface and downward in 
energy the Dirac cone (n-type) associated with the surface. This feature resembles the 
electronic structures found in graphene/ferroelectric/graphene composite [27], and gives rise 
to nontrivial spin textures for bands near the Fermi energy. For instance, the surface states 
of the P~ and P^ surfaces are expected to have the same spin helicity at the Fermi energy, 
as shown schematically in the inset of Fig.3b. This is confirmed by our first-principles 
calculations (Fig. 3c). Similar to non-ferroelectric topological insulators, currents on the P~ 
and surfaces carry opposite spins (Fig. 3d). Ferroelectricity allows these currents to be 
separated in energy as well. These results demonstrate that strained CsPbD is a FETI with 
TSSs strongly coupled to bulk polarization. 

Compared to a ferroelectric normal insulator (FENI), the EETI is expected to have more 
robust ferroelectricity. As shown in Fig. 4a, normal ferroelectrics (e.g., PbTiOs) are unstable 
under ideal open-circuit conditions (in vacuum) due to the large depolarization field. In 
practice, normal ferroelectrics are stabilized by electrodes or surface charges which eliminate 
the effect of the bound charges. An advantage of the FETI over the EENI is the guaranteed 
presence of metallic surface states, which can serve as innate metallic electrodes to create an 
intrinsic short-circuit condition, allowing for automatic surface charge dissipation and hence 
stable ferroelectricity in vacuum. 

The polarity of ferroelectric surfaces can be tuned and reversed by well-controlled ex¬ 
ternal stimuli, such as electric field, stress field and temperature. Therefore, the coupling 
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of switchable polarization and TSSs in FETI offers a powerful paradigm for novel device 
engineering. For example, the 180° domain wall in FETI may serve as a Veselago lens j2H] 
for precise focusing of electric current. Shown in Fig. 4b, the 180° domain wall separates 
domains with opposite polarities, thus creating an n-type Dirac cone on the surface and 
a p-type Dirac cone on the P~ surface. Perfect focusing of electric current occurs when the 
n-type and p-type Fermi circles have the same radius [29] . While this condition can only be 
obtained by careful tuning of local chemical potentials in graphene, it is automatically sat¬ 
isfied in FETIs because of the antisymmetric nature of the P^ and P~ surfaces. By moving 
the position of the domain wall via a gate voltage, it is possible to controllably change the 
focusing point on the P~ surface. 

In addition, ferroelectric materials can be patterned into designed multi-domain configu¬ 
rations. This enables simplified fabrication and control of high-density p-n junction arrays 
at the nanoscale (Fig. 4c), which may find applications in electronics and spintronics, such 
as electron beam supercollimation [HOI El] • The operating bandwidth of Dirac fermion-based 
electron beam supercollimators increases as the period of the p-n arrays decreases [SO]. This 
has proved to be a limiting constraint in designing graphene-based electron beam super- 
collimators using conventional metallic gates which are usually larger than 100 nm in size. 
In contrast, nanoscale ferroelectric domains can be easily created and manipulated [21E2], 
opening possibilities for spin-sensitive electronic devices of wide operating bandwidth. 

In conclusion, our ab initio calculations demonstrate that it is possible to realize ferroelec- 
tricity and robust topological surface states in a single material. The electronic structure 
of cubic CsPbIs with ferroelectric distortions at high pressure can act as a prototype for 
designing ferroelectric topological insulators that are stable at ambient conditions. Replac¬ 
ing Cs with larger cations (e.g., molecular cation, (CH 3 ) 4 N+) may help to suppress Pble 
octahedron rotation. The discovery of robust material realization of ferroelectric topological 
insulators can potentially lead to new device applications in electronics and spintronics. 


METHODS 

All DFT calculations are carried out with plane-wave density functional theory package 
QUANTUM-ESPRESSO [33] with Perdew-Burke-Ernzerhof (PBE) [3l] density functional 
and norm conserving, optimized [HS] , designed nonlocal [HE] pseudopotentials generated from 
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the OPIUM package. A plane-wave energy cutoff of 50 Ry is used. A five-atom unit cell 
is used for constructing the structural and topological phase diagram of cubic CsPbR with 
the Pblg octahedron rotation suppressed. An 8x8x8 Monkhorst-Pack fc-point grid is used 
to converge the charge density, while an 8 x 8 x 8 Monkhorst-Pack coarse grid, interpolated 
to a 32 X 32 X 32 fine grid, is used for calculations of the phonon entropy. A supercell slab 
method is used to calculate the band structures of surface states. Atoms in the slab model 
are fixed to their bulk atomic positions. The supercell for modeling the non-polar surface 
consists of 30 layers of ferroelectric CsPbR and 22 A of vacuum. The slab has surface normal 
along the [001] direction and the polarization along the [100] direction. The polar surfaces 
are modeled with a slab that has both surface normal and polarization oriented along the 
[100] direction. The vacuum inside the supercell is 15 A thick, and the dipole correction |37j 
is applied to remove the spurious interaction between different images. The surface band 
structure of polar surfaces is projected onto orthogonalized atomic wavefunctions, and the 
weight of atomic wavefunctions from P~^ and P~ surface atoms are then evaluated separately 
for 30 highest valence bands and 10 lowest conduction bands. 
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FIG. 1. I Pressure induced parelelectric—ferroelectric and trivial—non-trivial phase tran¬ 
sitions in CsPbla. a, Electronic band gaps of paraelectric and ferroelectric CsPbla as a function 
of lattice constant. Paraelectric CsPbIs has band gap closure at a = 6.19 A. Ferroelectric CsPbIs 
changes from a normal insulator (NI) to a topological insulator (TI) at a = 6.04 A. The system 
undergoes a paraelectric to ferroelectric structural phase transition as the lattice constant be¬ 
comes smaller than ps 5.90 A at room temperature, b, Evolution of band structure of ferroelectric 
CsPbIs as a function of lattice constant. At a = 6.04 A, a line node (not shown) appear near the 
time-reversal invariant point (i?) in the Brillouin zone, c, Electronic band gap and polarization 
for interpolated structures connecting paraelectric and ferroelectric CsPbIs with lattice constants 
fixed at 5.79 A. Along the path, structures have gradually increasing polarization and decreasing 
band gaps. The inset shows the adiabatic evolution of the band structure from paraelectric to fer¬ 
roelectric phase. The manifold of valence bands can be continuously transformed from paraelectric 
to ferroelectric without closing the band gap. 
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FIG. 2. I Electronic structure of non-polar surfaces of CsPbla. The surface normal is along 
the [001] direction. The shaded (grey) bands correspond to bulk bands. The red lines highlight 
the slab bands. The topological surface states emerge inside the bulk band gap. The tiny gap of 
the surface states (magnified in inset) appears due to the finite size of slab thickness 
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FIG. 3. I Electronic structure of polar surfaces of CsPbIs. a, Schematic of the slab geometry 
with surface normal and polarization along the [100] direction, creating P~ surface (left) and 
surface (right), b, Surface band structure. The red and blue lines represent the contributions of 
surface states from surface and P~ surface, respectively. The thickness of lines scales with the 
contribution. The inset illustrates the vertically shifted Dirac cones on the P^ and P~ surfaces, 
c, DFT spin texture for the lowest conduction band around M. The solid line denotes the F-M 
direction, d. Orientations of spin currents on polar surfaces. The spins on and P~ surfaces 
adopt opposite orientations for spin currents along the same direction. 
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FIG. 4. I Schematic of functional devices based on ferroelectric topological insulator, a, 

Robust ferroelectricity of FETI under open circuit condition. Ferroelectric thin film normal insula¬ 
tor (FENI) becomes paraelectric in vacuum due to the large depolarization field. The polarization 
of FETI is protected with metallic surface states serving as intrinsic metallic leads, b, Domain 
wall in FETI as Veselago lens. Electrons injected at S will focus at Di for domain wall positioned 
at line A. The domain wall can be moved to line B by reversing the direction of the polarization 
through a gate voltage. The focusing point will consequently be shifted to D2. c, Ferroelectric do¬ 
main arrays as nanoscale topological p-n junction arrays for electron beam supercollimation. The 
potential periodicity and amplitude can be tuned by varying the size of domain and the magnitude 
of the polarization respectively. 
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